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We examine the realization of a quantum CNOT gate by adiabatic operations. The principles of 
such systems and their analysis are briefly discussed and a model consisting of two weakly coupled 
double- potential well qubits is studied numerically. Regions of the parameter space with suitable 
well-defined sets of wavefunctions are found, in which then an adiabatic sweep of an external bias 
produces the switching behavior of CNOT. Results are presented on the adiabatic condition and 
the identification with the parameters of a flux-coupled two-SQUID system is given. For typical 
parameters adiabatic times in the nanosecond regime are obtained. 



The basic element of the "quantum computer" [1] is 
the quantum bit (qubit), a two level system, exhibiting 
quantum coherence between the states. Many physical 
realizations of the qubit have been proposed [2] . To ma- 
nipulate the qubit quantum gates [3] are necessary, logic 
devices capable of operating on linear combinations of 
input states. First there is the simple NOT, a one bit 
operation which can be viewed as an inversion operation 
on a qubit. The next step is to construct gates of a condi- 
tional character. A simple case to consider is the two-bit 
operation "controlled NOT" or CNOT. To realize such 
device it is natural to consider using an interaction be- 
tween the physical elements constituting the qubit. 

Among the possible mechanisms for manipulating cou- 
pled qubits adiabatic procedures are, as explained below, 
of special interest. Furthermore, it has been suggested 
that adiabatic procedures may be robust with respect 
to certain kinds of errors [4]. In particular with super- 
conducting devices, Averin [5] has suggested using small 
Josephson junctions in the coulomb blockade regime, and 
we have mentioned the possibility of using SQUID qubits 
with flux coupling [6]. In this letter we will explain 
some general principles for studying such systems and to 
present numerical calculations relevant to their behavior 
and design. 

CNOT is a two-qubit operation and we will represent 
it by two interacting double-potential well systems. Each 
double well system may be though of as an approximately 
independent qubit since we shall keep the coupling weak. 
Qualitatively, we will use the procedure of performing 



an adiabatic NOT [6] on the first qubit while trying to 
influence its behavior by the state of the second. We find 
a region of parameter space where this works. 

Hamiltonian: We take the following model hamiltonian 
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containing dimensionless 
"position" coordinates <j>\. 
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4>2 , while the potential term is 
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The function f(4>) is chosen so that a double-well po- 
tential results for each <j> variable. We shall use f(d>) = 
1 — \<f> 2 + ^(j) 4 '-. Vo, h, h, h2, Pi and /3 2 are constants 
depending on system parameters. The two 4> ext are exter- 
nal biases which can be adjusted or varied to perform the 
operation and to find favorable operating points for the 
device. Their values determine the degree of asymmetry 
of each double well system; when they are zero the wells 
are symmetric (for 1x2 = 0). Note that when the cou- 
pling parameter i 12 is zero the hamiltonian simply rep- 
resents two non-interacting systems. Our model system 
thus consists of two weakly interacting double-potential 
wells with externally adjustable biases 4>l xt , <l>2 Xt - Fig 1 
shows the equipotential contours of V((f>i, $2), with its 
four potential wells. 

Representation of Logical States: We first require a 
representation of the four states of the two-qubit system. 
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FIG. 1. Potential as 
in Eq [2], with its four 
wells. The coordinate 
4>\ (target bit) runs hor- 
izontally and fa (control 
bit) vertically. Red indi- 
cates the deepest poten- 
tial, purple the highest. 



Each logical state will be represented by a wavefunction 
localized in a distinct potential well. When this obtains, 
the logical state corresponds to a distinguishable physi- 
cal state with high probability. With SQUID qubits, for 
example, if for a single SQUID a wavefunction concen- 
trated on the left of the double potential well corresponds 
to current going clockwise, then for a two-SQUID system 
the lower left well of Fig 1 corresponds to current counter- 
clockwise in SQUID 2 and current clockwise in SQUID 
1. Since, as will be explained in the next paragraph, we 
work with energy eigenstates, a first requirement on the 
hamiltonian is thus that it yield a set of "good wavefunc- 
tions" , that is where the first four energy eigenstates are 
well localized in the four different wells of Fig 1. 

Assuming such a set of wavefunctions has been found, 
visualization of the situation is aided by the use of 
tableaux indicating where the wavefunctions are local- 
ized on Fig 1. Labeling the first four energy eigenstates 
in order of increasing energy 1,2,3,4, examples of possible 
tableaux are seen in Eq[3]. 

Representation of CNOT: Our aim is to represent a 
logical operation such as CNOT by a mapping of the 
set of initial logical states to a certain set of final logical 
states. This will be represented by a particular rearrange- 
ment of states on the tableau. 

CNOT is defined by the conditions: A) the control bit 
does not change its state, and B) the target bit is reversed 
or not reversed, according to whether the control bit is 
1 or 0. If we identify the top row of the tableau with 
control bit =0 and the bottom row with control bit = 1, 
a physical embodiment of CNOT would be 



4 3 
1 2 



4 3 
2 1 
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Condition A) on the stability of the control bit is exhib- 
ited in that no states move between the top and bottom 
row. Condition B) is realized in that the top row remains 
unchanged while the bottom row is "flipped" . 

Adiabatic Operations: Realization of operations such 
as Eq[3] can be accomplished in an especially transpar- 
ent way by using adiabatic processes. This is due to the 
"no level crossing" behavior of adiabatic evolution. The 
no-crossing property assures that a state initially in the 
first, or second, or third,.... energy level will end up in 



the first, or second, or third,... energy level after the 
adiabatic evolution, while at the same time the physi- 
cal properties associated with the level may be changing. 
Thus in Eq[3] with SQUIDs, state 1 begins as a config- 
uration with the current clockwise in SQUID 1, counter- 
clockwise in SQUID 2 and ends up as a configuration 
where the current remains counter-clockwise in SQUID 
2 but is now reversed to counter-clockwise in SQUID 1. 

One can proceed as follows: we search for an initial 
hamiltonian whose variable parameters (<j>l xt , 4>2 Xt ) are 
adjusted to give the left tableau of Eq[3]. Then, we search 
for a final hamiltonian where another set of {<fii xt ,<f>2 Xt )i 
gives the tableau on the right. If the two parameter sets 
can be connected by a smooth, slow transformation,-a 
"sweep"- we have obtained an adiabatic realization of 
our operation, here CNOT. 

In this procedure we need only to study the station- 
ary Schroedinger equation at first. This is an important 
simplification for the numerical analysis. However, after 
having determined some suitable parameter sets we shall 
also study the full time-dependent Schroedinger equa- 
tion. This is necessary to determine what sweep speed is 
"slow" , that is guarantees adiabatic behavior. 

Numerical Methods: Our problem involves two vari- 
ables and tunneling through four barriers, as well as 
a multi-dimensional parameter space. To deal with 
this complex situation we turn to a recently developed 
method [7] for numerical solution of the Schroedinger 
equation. A large basis of harmonic oscillator wavefunc- 
tions is used to reduce the problem to an array of fast 
algebraic manipulations, programmed in Mathematica. 
Except for the small Z12 ~ 10~ 3 , we work with parame- 
ters of order one, hence the resulting dimensionlcss ener- 
gies are also of order one. However the splittings among 
the lowest levels, which are what we manipulate, result 

from tunneling and are small - 10 -3 1(T 4 . Hence 

four place accuracy is necessary. Using these methods we 
have been able to find a region of the (<fil xt , 4>2 Xt ) param- 
eters space where there are "good wavefunctions" . These 
are indicated as the gray regions of Fig 2, with the pa- 
rameters as indicated. Reducing the value of Z12 leads to 
a shrinking of these regions on the plot. 

As a by-product of our numerical work we can also 
examine the validity of the frequently used "psuedo-spin" 
picture. One often usefully visualizes [6] the lowest quasi- 
degenerate levels of the system as "spins" . This picture 
requires, however, that the moving statevectors remain 
in the hilbert space spanned by an initial (here four ) set 
of states. By evaluating wavefunction overlaps we find 
this is true, to a good approximation, supporting the use 
of the "spin" picture. We stress, however, that we do not 
need this simplification in our calculations. 

Switching Behavior: We have been able to obtain 
switching behavior according to Eq[3] for the "good" re- 
gions of Fig 2, by means of the following operation: the 
control bias <j>2 Xt is held constant at a relatively high value 
while there is a sweep of the target bias 4>l xt . This is 
a generalization of a simple NOT [6] on <f>\. The results 
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may be understood in terms of a simple model: the I12 
coupling produces an extra bias on the target bit which 
"helps or hinders" the NOT operation. 

The relatively large bias on <j>2 comes from condition 
A): we attempt to "immobilize" the control bit despite 
the perturbations communicated by the sweep of fa via 
l\2- We therefore investigate the region \fa xt \ << l^'l- 
If (j>2 is indeed successfully "immobilized" it will be fixed 
in one of its two potential wells and can have only the 
values <j>2 ~ ±1. As seen by fa, this amounts to ± an 
extra bias. To linear order (since we take all fa xt , fa xt 
small compared to 1 and fa , fa are in the neighborhood 
of 1) and introducing the notation 

4>n ff =<t>? t ± 1 -^ (4) 

the potential terms involving fa in Eq [2] become 

-2hfafa Xt -2l 12 fa(±l) 

= -21^ ± l f) = -ihwrtff . (5) 

As a consequence, there is an effective shift in the ex- 
ternal bias on fa by (±7^ )■ This is just as suggested 
by the "help or hindering" picture, and we further learn 
that the magnitude of the "help" is j^-. According to 
this picture we should try to analyze the behavior of fa 
(target bit) as if it were simply under a modified external 
bias 4>f X eff- 

We would then expect if a certain tableau obtains and 
fa xt is varied, that the tableau is maintained until there 
is a sign switch for fa x lfj- To see how the sign switch 
occurs, note that for fa xt the sign of ^Tlff IS given 
by the ± from the control bit. But for \fa xt \ » 
the sign is controlled by fa xt itself. There is, therefore, a 
region around fa xt = where the tableau is determined 
by the control bit, and another region for large \fa | 
where the tableau is determined by the sign of fa xt . As 
we cross from one region to another, one pair of states will 
retain the sign it had for large \fa xt \ and the other pair 
of states will switch. This would be the desired behavior. 
According to Eq [4] the point where the switch from one 
tableau to another should take place is given by \ fa xt \ ~ 
I12/I1. 

Numerical Results: The three different gray areas of 
Fig 2 have well defined tableaux as follows: 

(? i) 5) (2 ;) ■ («) 

for the intermediate gray region (left), the darkest region 
(center) and the light gray region (right), respectively. 
These tableaux fit with the description arrived at in the 
model, where either the top or bottom row inverts as we 
go from the central region to large l^f^l. Hence a sweep 
from the central region to the right region produces the 
mapping Eq [3] . Similarly sweeping from the right region 



to the center and from the left region to the central re- 
gion and vice- versa can also serve as realizations, differing 
simply in the assignment of (0,1) for the quantum states 
or the names for the qubits. Fig 2 shows that the switch 
between tableaux occurs quite close to \fa xt \ = \li2/h\, 
as predicted by the model. Also, it seems the inequality 
\fa xt \ « \fa xt \ need not be very strong for "immobiliza- 
tion" . 

Adiabatic Condition: An important time scale is Tadiab-, 
the shortest time in which an operation can be performed 
adiabatically. This time is relevant both to the maximum 
speed of the device and with respect to decoherence and 
relaxation effects since the operation must take place in 
times short compared to those for these effects. We thus 
now examine the time dependent Schroedinger equation 
iJ^ip = Hip. We work with the dimensionless time vari- 
able r, where the connection to usual time t is given by 
t = Eot. Eq is an energy parameter in electron volts or 
Hz, characteristic of the particular system under consid- 
eration. It also gives the overall energy scale. Thus all 
energies are in units of Eq and all times are in units of 

Eo 1 (ft = I)- 

An estimate for NOT [6] gave 

Tadiab = euJ tunnel ~ eT rabi i CO 

where e is the length of the sweep in energy and u^~ unne t = 
T ra i,i , the inverse tunneling energy or oscillation time be- 
tween the two states at minimum separation. Since here 
we also perform a kind of NOT, we expect a smiliar 
relation to hold, where u>tunnel or T~ abi is the smallest 
level splitting during the adiabatic passage and the sweep 
length e may be read off from the total energy shift of 
the wells. We define a "degree of adiabaticity" by per- 
forming a sweep numerically and taking the overlap of 
the resulting wavefunction with the wavefunction of the 
corresponding stationary final hamiltonian; that is to say 
the overlap with the wavefunction that would result from 
an infinitely slow sweep. The square of this amplitude, 
Pfi, is shown in Fig 3 as a function of sweep time, for a 
sweep (0.002,-0.01) -> (0.008, -0.01) on Fig 2. The ar- 
row indicates the theoretical prediction using Eq[7]. As 
would be expected, r a di a b ~ e 6 2 w 500 is a large number 
in dimensionless units. 

Identification with SQUID parameters: Eqns [1,2] arise 
by standard methods in the analysis of two rf SQUID 
loops coupled by a mutual inductance L12 [8]. The 
Josephson relation leads to f(fa) — cos0, to which our 
/(</>) = 1 — \<t> 2 + ^j'/' 4 is a good approximation. 

With L\, L 2 the SQUID inductances and C\, C2 the 
Josephson junction capacitances, the energy scale factor 
is Eq = 1/s/LC, where C = and L = Ll ^~ L ^ » 

yJL\L2. A set of reasonable values for the SQUID is 
Li = 300pH, L 2 = 280pH, L 12 = 1.8pH, d = C 2 = 
O.lpF and (3\ = (3 2 = 1-28. Since in frequency units 
one finds En w . 1 1000 GHZ these values give 

^/L/pH C/pF 5 

Eq w 185 GHZ. Then the T a di a b ~ 500 in dimensionless 
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FIG. 2. A region of the <j>\ xt , cj>l xt plane with a well de- 
fined set of wavefunctions as explained in the text. The 
coupling parameter is Z12 = 0.005. Other parameters were 
h=h = 1, Pi=fo = 1.19, fn = fi 2 = V = 16.3 
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FIG. 3. The adiabaticity parameter versus sweep time for a 
sweep of (4>\ xt , 4>% xt ): (0.002, -0.01) -> (0.008, -0.01). P fi =l 
denotes perfect adiabaticity. The arrow indicates the theoret- 
ical estimate T adiab = eui^ nnel . 
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units corresponds to T ac iiab ~ 500/E rts 2.7 • 10 9 s in 
seconds. 

Finally we note that here, as with all discussions of 
quantum computation, the important open question is 
the time scale for decoherence Tdec- In [6] we estimated 
T dec « {few) x 10 _6 s for the SQUID at 40mK (and also 
suggested a method for its direct measurement). With 
the above estimate for T a diab it thus appears, at least for 
the SQUID, that the adiabatic condition allows for opera- 
tions in times less than Tdec- However, the question of the 
decoherence time is controversial and system-dependent 
and will probably only be resolved convincingly by ex- 
periment. In this respect, it is encouraging that in ex- 
periments [9] showing evidence for macroscopic quantum 
behaviour of the SQUID the dissipation value is rather 
small. 
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